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Abstract 

The QCD+QED axial anomaly low-energy theorems are applied to estimate the 
accuracy of an instanton vacuum - based chiral quark model. The low-energy the- 
orems give an exact relation between the matrix elements of the gluon and photon 
parts of the QCD+QED axial anomaly operator equation. The matrix elements be- 
tween vacuum and two photon states and between vacuum and two gluon states are 
calculated for arbitrary Nf. It is shown that this model does satisfy the low-energy 
theorems with an accuracy of ~ 17%. 

We estimate also the contribution of the nonperturbative conversion of gluons 
into photons to the decay rf — > 27 and compare with experimental data. 

PACS number(s): 11.15.Tk , 12.38.Lg. 
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1 Introduction 



In gauge theories the symmetries of the initial classical lagrangian may be destroyed 
by quantum fluctuations. One of the most important examples is the famous axial 
anomaly. The axial anomaly leads to many interesting nonperturbative phenomena 
in physics. Among them are S-violation processes in electroweak (EW) physics, 
Ua(1) problem in QCD etc. The solution of these problems is intimately related to 
the topologically nontrivial structure of the vacuum in the gauge theories. 

In this paper we apply the axial anomaly low-energy theorems to test the chiral 
quark model which is based on the instanton model of QCD vacuum. We find that 
this model does satisfy these theorems with an accuracy of ~ 17%. This conclu- 
sion provides solid background to calculate different amplitudes of nonperturbative 
conversion of gluons into hadrons and photons. 

We estimate also the contribution of the nonperturbative conversion of gluons 
into photons to the decay rf —>■ 27 . By using experimental date for the width 
of rf — > 27 we conclude that //-singlet axial current coupling constant is given as 
/o 

1.1 Low-energy theorems from axial anomaly 

The axial anomaly in the divergence of the singlet axial current in QCD + QED 
leads to a low-energy theorem for the matrix elements of this operator equation over 
vacuum and two-photon states: 

(0|jv >^ g6|27) = N ^Y,Q 2 f Fil)p{2) C 1 ) 

at q 2 = H] . Here, Nf is the number of the flavors, g is the QCD coupling constant 
with GG = -e^"^ T G a G\ a and G a being the operator of the gluon field strengths , 
iV c is the number of the colors, e and Q f are QED coupling constant and the electric 
charges of the quarks respectively, 2F^F^ = e^ Xf7 F$F ( *J, F$ = efq iv - e®q ifi , 
e^' 2 \ qi t 2 are polarizations and momenta of photons respectively and q = qi + q2- 
This relation is a consequence of the absence of a massless singlet pseudoscalar 
boson. Here the contribution of the quark masses is neglected. 

As it is evident, gluons can interact with photons only through quark loops. In 
perturbation theory it leads to at least ~ g 4 results for the left hand side of Eq. ([!]) 
(see e.g. recent discussion of the higher- loop contributions to the axial anomaly M). 
So, the solution of this theorem is related only to the nonperturbative phenomena 
connected with the structure of QCD vacuum. 

Another nontrivial low-energy theorem concerns the matrix element over vacuum 
and two-gluon states: 

(0\g 2 GG\2gluons) = (2) 

at the limit q 2 = 0. 

These matrix elements are calculated in the instanton vacuum generated N-JL 
type quark model J|, ||J] for arbitrary Nf. 

1.2 The instanton vacuum of QCD 

The instanton is the solution of gluodynamics in the Euclidian space[jS|: 
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A%{x) - 2g X Off\l v[{x _ P z)2 X + ^ x _ z)2 . (3) 

The anti-instanton solution has the same form as in Eq. (|]) but with rf instead of 
fj^ v . Here is the orientation matrix of the instanton I (anti-instanton I) in color 
space, 77(77) - t'Hooft factors 0, p is the size and z is the position of the instanton. 
For large interinstanton distances R» p the sum of the N + instantons and the N_ 
anti-instantons is also an approximate solution. The calculation of the action for 
gluon fields leads to a sum of actions of free instantons and classical interinstanton 
potential V(R, pi, P2, O), where O = OJO2 is a matrix of relative orientation. The 
most important part is the instanton-antiinstanton potential Vjj. It is well-known 
that at large distances it has a form of the dipole-dipole interaction potential and 
may be attractive. 

The main assumption of the model concerns small distances R ~ p. At small dis- 
tances it is assumed that there is a repulsion. This assumption is supported by both 
phenomenological and theoretical considerations 0,0. This leads to a stabilization 
of the size and density of instantons. The quantities mean number of instantons (N) 
and a mean instanton size p include effects of the instanton interactions Q. In the 
following calculations, for simplicity, all instanton sizes are considered to be p. 

Both the phenomenological estimates and variational calculations lead to a mean 

, x 1/4 

interinstanton distance of R = yjj^J ~ 1 fm and a mean instanton size of p ~ 
1/3 fm. The small packing parameter (p/R) 4 = 0.012 provides a possibility for the 
independent averaging over positions and orientations of instantons. 



1.3 The chiral quark model 

The main assumption of this model is an interpolation formula for the quark propa- 
gator in the single instanton field. It is approximated as the sum of a free propagator 
and an explicit contribution of the zero mode ||, 

UV(£ J( i))+OTi) . «0<9) - ^ 4 

\ v ' / 1— inst ™, 



Here, V == d — igA, <&±(x; £im) is the zero mode wave function of the fermion in 
the background of one J(/)|J. It depends on the collective instanton variables £jrj\ 
- the size p, the position z and the orientation O of the instanton. 

This interpolating formula should be accurate both at small momenta (p l/p), 
where the zero mode is dominant, and at large momenta (p ^> l/p), where the prop- 
agator reduces to the free one. In the background of an iV±-instanton configuration 
and keeping in mind the low density of the instanton media this formula leads to 
the partition function of the model: 



J DipDvp^ exp( J d^xipHdip) W+ + W*~, (5) 



where 



J ± (z) f9 = j e M-*(k ~ l)z) F{k)F{l) 4(k)i(l± 75 )^(/), (7) 
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and the contribution of the current quark masses is neglected. The form-factor F 
is related to the zero-mode wave function in momentum space $±(fc; £/(/)) and is 
equal to: 

F(k) = -tj f [I (t)Ko(t) - h(t)Ki(t)] - { l r *3^ °_ ^ , (8) 



with t = \kp. 
The formula: 



(ab) N = JdX exp{Nln^- - N + Xb) (9) 
(here N >> 1) provides the final expression for the partition function ||: 

Z N = J DtpD^exp( J ijHdtp + Y + + YJ), (10) 

where 

r /2V\ N f~ 1 r 

Y± = (i) N f\ J d 4 z det J±{z) = ( — J (iMff J d 4 z det J±(z). (11) 

The self-consistency condition at the saddle point in Eqs. (]10|) , and fllT] ) leads to 

r d A k M 2 FHk) , , 

an V / — — = N (12) 

c J (27r) 4 M 2 F 4 (fc) + F { ' 

2 Calculations with the low-energy theorems 

In the quasiclassical(saddle-point) approximation any gluon operator receives its 
main contribution from instanton background. As an example, for one instanton 
(anti- instanton) I (I), 



9 2 G\x)= 192P ^Ji.r- :) (13, 

[p 2 + {x — zY\ 

and 



g 2 GG(x) = ±f(x-z). (14) 

In this particular case, the calculation of the matrix element can be reduced to 
the calculation of the partition function 

Z n [k,o\ = Zn 1 J DipD^ exp(-S ef f), (15) 

with an effective action, S e ff, in the presence of an external electromagnetic field, 
a^, and an external field k(x) is given as 

-S eff = fipHDiP + Y + + Y_ + fdx (Y GG+ (x) + Y GG _(x)) k(x), 

(16) 



Nf — l 



X) K(X) 



Y GG± (x) = ± (f ) 1 (iM) N f J d 4 z f(x — z) det J±(z), 
where D = d — ieQfd. Finally, Eq. ( |i~6"|) can be rewritten as: 
-S eff = j^W^+ 

(17) 

(¥ Y^fdz det(tM + (z)J + (z)) + (f) ' 'jdz det(iM_(z)J_(z)), 
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where 

M±{z) = f 1 ± / dxn(x)f(x -z)\ M. 
Another remarkable formula 

exp(Adet[L4]) = J dMexp -(N f - l)A _w F T (detA4) w F T + itr(MA) 

is used in the following discussion. It is possible to check this by the saddle point 
approximation of the integral. 

By using Eq. (|T8D it is easy to show that the effective bosonized action S e ff[A4±, a, k), 
describing mesons (the matrices M.±) in the presence of the external fields a M and 
k is 

2V X ~ 1 



-S eff [M±,a, K ] = Jdz {-(N f - 1) (f j (det M ± ) N f-' J + 

Tr In [if) + iM+MF 2 (1 + (nf)) N f * |(1 + 7s ) + iM-MF 2 (1 - («/)) V |(1 - 75 ! 

(19) 

For the processes without mesons(A^± = 1) the partition function is: 

Z N [K,a]=ex P Trln(tb + tMF 2 (l + (Kf)) N 7 1 ^(l + l5 ) + 
iMF 2 (l-( K /))Vi(l- 75 )) x (id + iMF 2 )' 1 , 

where («/) = / dxn(x)f(x — z). 

2.1 The low-energy theorems for the matrix element be- 
tween vacuum and two- photons states 

The matrix element is generated by 

8Zn[k, a] 



Sn(x)5a ll (xi)5a u (x 2 ) 



|ft,a=0 - 



As it is clear from Eq. ( p£|) the factors in the vertices in the corresponding Feynman 
loop-diagram are and iMfF 2 'j 5 Nj' 1 . 

We must calculate A(g 2 ) (in the limg 2 — > 0), which is defined by: 

(27i) 4 5(q - qi - q 2 )A(q 2 ) = J dx exp(-iqx)(0\g 2 GG(x)\2 1 ) = 

e 2 eW4 2 ) J(0\T(g 2 GG(x)jf l m (x 1 )jt m (x 2 )\0) expi(-iqx + q x x x + q 2 x 2 )dxdx x dx 2 . 

(21) 

It is clear that 

A(q 2 )=^^f(q 2 )N c e 2 j:Q 2 f 

f 

r d 4 p »MF!F 2 7 5 (p - g t + tMF^T^p + iMF 3 2 ) 7t; (p + g 2 + »MF 2 2 ) 
X n J (2vr) 4 ((p- qi ) 2 + M 2 Ff)(p 2 + M 2 Fi)((p + q 2 ) 2 + M 2 Fi) 

+ ( fi ^u,q 1 ^q 2 )). (22) 
Here F x = F((p - F 2 = F(p + q 2 ), F 3 = F(p) and 

f(q 2 ) = J dxexp(-iqx)f(x) (23) 



Table 1. The accuracy of the model for the different values of the parameter 
r = Mp in the Eq. fl2"7p to satisfy axial anomaly low-energy theorems. 



r 


.0 


.1 


.2 


.3 


.4 


.5 


.6 


.7 .8 


.9 


1.0 


5{r) 


.00 


.05 


.09 


.13 


.15 


.17 


.19 


.20 .21 


.22 


.23 



is the form-factor of the one-instanton contribution to g 2 GG. At q 2 = 

f(q 2 = 0) = 327T 2 . 
It is easy to show that the trace in Eq.(p2|) can be reduced to 

8M 2 e^q lx q 2a T(q 2 ,r), 

where 

IV r) = f ^3 + (24) 

W ' J 7 (2tt 4 ) ((p - gi ) 2 + M 2 F 4 ) {p 2 + M 2 F 4 ) ((p + q 2 ) 2 + M 2 F 4 ) 1 ' 

and r = Mp. It is possible to calculate this integral analytically if we put r = in 
Eq. (|24D (In this case F — 1). In this approximation and at a small values of q 2 , we 
have 

r ^ r = °» = «( 1 -i^) (25) 

As a result, in the r = approximation the left side of the low-energy theorem, 
Eq.®, is 

( N f)(^iYQi) F m F w (26) 

and this coincides with the right side of Eq. ([I]) . 

The calculations of the integrals in Eq. ( p4| ) can be performed also with an 
account of the form-factor F(p). The deviation of results of the calculation from 
the approximate results of the r = tell us about the accuracy of this model. We 
need to calculate 

rtf = o, r) = tf //>«^g) (27) 

V 7 y (2tt 4 ) (p 2 + M 2 F 4 (p)) 3 V ; 

An evolution of T(q 2 = 0, r) was performed numerically. Table 1 presents the factor 
5(r) = 1 - T(0, r)/T(0, r = 0), where r = Mp. 

2.2 The low-energy theorem for the matrix element be- 
tween vacuum and two-gluons states 

Here we present, without details, calculations related to Eq. (§). This matrix 
element can be written in the form: 



\g 2 GG\g(eW, qi ),g(eV,q 2 )) = e«^eg^ 



(28) 



x Jd 2 d 2 $\Tg 2 GGA%{xx)A%(x 2 )\Q) expi(q lXl + q 2 x 2 )dx 1 dx 2 . 
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Here A a ^(x) is a total gluon field, e^ a \ qi are the polarization and the momentum of 
gluons respectively. 

As usual, we expand the total field A^(x) around the instanton background. 
The main term in Eq. (|28|) is the contribution of the instanton background and is 
~ 0(g~ 2 ). The next term is the contribution of the perturbative fluctuations over 
instanton background and is ~ 0(g 2 ). It is easy to see from previous considerations 
that 0(g~ 2 ) term is given by the formula 

Z N l J DipD-ip^ eMSeff) ((Y GdAA+ (x) + Y GGAA _(x)) Q) , (29) 

where 

Yggaa± = ± (f (iM) N f J d±z f{x - z) 

(30) 

x JdO(-d 2 )A^(x 1 )(-d 2 )A I ^(x 2 ) det J±(z), 

Here the instanton(anti-instanton) is located at the point z with its orientation O. 

Repeating the bosonization trick leads to the result for the 0(g~ 2 ) contribution 
which is proportional to 

Tr[(id + iMF 2 )~ 1 iMF 2 >y 5 \. 

It is clear that this Tr and as a consequence 0(g~ 2 ) term are equal to zero. 

The next order (0(g 2 )) term is the contribution of the two diagrams. The first 
diagram is the direct contribution of the operator g 2 GG which is equal to 

-g 2 G^& 2 \ 

where 2G^G^ = e^G^ a G^ a , G$ a = efq lu - e®q i(l . 

The factors in the vertices of the second loop-diagram are gX a /2 , ~f^ and iMfF 2/ j 5 N 
An account of the contribution from all flavors gives the coefficient Nf. 

A comparison with the previous calculations (Eq. (|2ll) , (p2"l)) leads to the result 
that the contribution of the second loop-diagram is equal in magnitude but opposite 
in sign to the contribution of the first diagram at q 2 = 0. 

Then, terms of 0(g 2 ) are equal to zero in the limit q 2 — > 0. This conclusion 
is valid in r = approximation. Thus we conclude that the instanton vacuum 
generated chiral quark model satisfies the low-energy theorems with an accuracy 
given in Table 1. 

3 Nonperurbative conversion of gluons into pho- 
tons and rj' —> 27 decay 

One of the first applications of QED axial anomaly was electromagnetic decays of 
pseudoscalar mesons |J. Here we consider the application of the QCD+QED axial 
anomaly to the same kind of decays. Let us consider the matrix element of the 
divergence of the axial singlet current <9 M j^ between vacuum and two-photons states 
at nonzero but small q 2 . We will compare this one with similar matrix element of 
the divergence of the third component of the axial isovector current d^j^ 3 . The 
definition of currents are j b = q^^q and j^ 1 = ^7^75(^/2)5. The matrix elements 
of these currents < 0\j^\r]'(p) >= ^^ifop^ and < ^\j^ l \^j{p) >= Sijif-^p^ defines 
the couplings fo and 
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The decay constant of charged pions is = 93.3MeV. In contrast to this one, 
fo has not been measured experimentally. Our approach provides a method for 
estimating f from the data on r/' — > 27 decay 

The matrix element under consideration is 

(0|^|2 7 > = ( |JV>^G?G|27> -N C £- 2 ^Q}F^F^ (31) 

We neglect here the masses of the u, d, s quarks. 

The lowest intermediate meson state in this equation is 7/ - meson with the 
non-zero mass m v / due to axial anomaly, as was mentioned above. With definition 
(r/'l 27) = fn'-^i^F^ 1 ' it is easy to conclude that 0(q 2 ) terms provide the equation 

In the r = approximation we conclude that right side of this equation is equal to 

^^(ji^A (33) 

On the other hand, as it is well known, the QED axial anomaly in the divergence of 
the third component of the isovector axial current gives the expression 

Uf^ = N c ^(Ql-Q 2 d ), (34) 

where fn ->2-y is defined analogously with Theoretical value of the width 

^0^27 = 7.25eV is rather close to the experimental one 7.95eV. 

With experimental value for the width T v ^ 2l = 4:.296KeV, M = MOMeV and 
p = 0.3/m in fl3"3] ) we estimate that 

fo = 1.16/ w . (35) 

This quantity was estimated previously by another method using 5'?7/(3)-symmetry 
type arguments and assumptions about the J/ip — > 7*7(7/') decay mechanism [|K| 
[TTj, [T3]|. Here the axial anomaly equation, eq. (|3"TD, is used again but for the 
matrix elements over the vacuum and 77 and 77' states. As a result, neglecting the 
contributions of the current masses of the light (u, d, s) quarks we get 

m^fo = N f ^- 2 < 0\GG\rf > (36) 

The matrix element on the right side can be extracted from data on the transi- 
tions and decays of heavy quarkonium systems fTT|j . Taking into account theoretical 
relation [JR], [TT| and experimental widths of the J/ip — > 777(7/') decaysj 

I < 0\GG\ V > > \ 2 P 3 V , _ r(j/v> - rfM) _ - 



I < 0\GG\rj > | 2 V% r(J/V -> 77/(7),)) 
it is easy to find that |10| 

1 — M — = 2.46. 

< 0\GG\r] > 
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On the other hand 

N fl7Tl < °\ Gd ^ >= - 2im ° < °\^s\v >= (h l/2 fsm 2 v . 
L6n z 3 1 

577/(3) symmetry arguments give f$ = 1.86 f n . 
This type of the estimates gives 

fo = l-06f n . (37) 

The estimates, Eqs.(p5|) and (|37 ) agree with each other even better than the usual 
accuracy of the S77f(3)-symmetry. It is crucial to repeat all of these estimates with 
account of the 5£7/(3)-symmetry breaking. 



4 Conclusion 

The solution of the axial anomaly low-energy theorems very nontrivially related 
with nonperturbative instanton structure of the QCD vacuum. The accuracy of the 
instanton vacuum based chiral quark model || to satisfy these theorems is ~ 17% 
and, probably, is related with the accuracy of the interpolating formula (f|) for the 
quark propagator. 

This approach provides a solid background to calculating the different amplitudes 
of the nonperturbative conversion of gluons into hadrons and photons. As example, 
we applied this approach to the rj' — > 27 decay process and, in result, successfully 
estimated ^'-singlet axial current coupling constant, the Eq.(|35|). 

We are planning to apply this approach to the hadronic and photonic transitions 
in a heavy quark systems and the scattering of these systems on hadrons. 

We hope to provide further discussions of all of these problems in a separate 
publication. 
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